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ABSTRACT 

Planets form in the discs of gas and dust that surround young stars. It is not known whether gas giant 
planets on wide orbits form the same way as Jupiter or by fragmentation of gravitationally unstable 
discs. Here we show that a giant planet, which has formed in the outer regions of a protostellar disc, 
initially migrates fast towards the central star (migration timescale ^ 10 4 yr) while accreting gas from 
the disc. However, in contrast with previous studies, we find that the planet eventually opens up a 
gap in the disc and the migration is essentially halted. At the same time, accretion-powered radiative 
feedback from the planet, significantly limits its mass growth, keeping it within the planetary mass 
regime (i.e. below the deuterium burning limit) at least for the initial stages of disc evolution. Giant 
planets may therefore be able to survive on wide orbits despite their initial fast inward migration, 
shaping the environment in which terrestrial planets that may harbour life form. 

Subject headings: protoplanetary disks - planet-disk interactions - planets and satellites: formation, 
gaseous planets - hydrodynamics 


1. INTRODUCTION 

Observations of exoplanetary systems are important 
as they provide significant clues about the formation of 
our Solar system. In the last few years direct imaging of 


massive giant i 

planets on wide orbits has become possible 

(iMarois et al.l 

20081; [Ireland et all 120111: 

Kuzuhara et al. 

2013 

: lAller et al.ll2013: Bailev et al. 20U 

t iRameau et al. 

2013 

; Galicher et al. 2014; Kraus et al. 

20141). These 


planets challenge our understanding of planet forma¬ 
tion as they are unlikely to form by core accretion, 
i.e. the coagulation of dust particles to progressively 
larger a ggregates and an accretio n of a gaseous en- 
velope (Safronov & Zviagina 1969; Goldreich & Ward 
Il973t iMizunol 1 19801: iPollack et al.l Il996h . Gravitational 
fragmentation o f protostellar discs is an alternative for¬ 
mation scen ario (jKuiperll 19511 : lCameronlll978l: lBosslll997t 
lBolevll2009l ). It has been argued that protoplanets form¬ 
ing in gravitationally unstable discs quickly migrate to¬ 
ward s the host star (jBaruteau et all 120 111 : iMichael et al.l 
12011 . where they may be disrupted by tidal torques 
([Navakshin fe Cliall2Q13t iTsukamoto et al.ll2013l) . More¬ 
over, it has been suggested that protoplanets rapidly 
grow in mass by accreting material from the disc to be¬ 
come brown dwarfs and low-mass stars ([ Rafikov 2005; 
iStamatellos fe Whitworth!l2009at IZhu et al.ll2012f ). 

In this Letter we present models of the interaction of 
a newly formed massive planetary embryo (henceforth 
protoplanet ) with a protostellar disc that contradict the 
above claims. We critically include in our models two in¬ 
gredients that their combined effect has been ignored by 
previous studies: (i) gas accretion onto the protoplanet, 
and (ii) accretion-powered radiative feedback from the 
protoplanet. 

2. COMPUTATIONAL MODEL 

We assume that the protoplanet has formed in the pro¬ 
tostellar disc of a Sun-like star by gravitational fragmen- 
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tation. For this to happen the disc has to be relatively 
massive and the protoplanet’s initial orbital radius is 
at least 50 — 100 AU, as the conditions are favour able 
for fragmentation only in the dis c outer region (jRafikovI 
120051: IStamatellos fe Whitworth] 12009a ). The mininum 
initial mass of the protoplanet is set by the opacity limit 
for fragmentation which is though t to be 1_—_5 Mj, 
where Mj is the mass of Jupiter dLow fe Lvnden-Belll 
19761: 1 Whi tworth fe Stamatellosl 120061: iBolev et al 
Kratter et al nHHI. Here we examine the evolution of 
such a protoplanet as it interacts with its parent disc, 
using a 3D Smoothed Particle Hydrodynamics method 
that includes the effects of radiative transport in the 
disc (|Stamatellos et al.l 1200 7: IStamatellos fe Whit wort hi 
I2009bl) . 

2.1. Initial conditions 

We assume a star-disc system in which the central star 
has initial mass M* = 1M 0 . The initial disc mass is 
M d =0.1 M 0 and the initial disc radius is R n = 100 AU. 
The disc is represented by 10 6 SPH particles. The inner 
disc boundary is set at 0.2 AU. The disc initial surface 
density is 

E 0 (iZ) = E(lAU) (dU) , (1) 

and the disc temperature 

/ d \ —3/4 

r o (i?)= 250K (—J + 10 K, (2) 

where E(1AU) is determined by the disc mass and radius, 
and R is the distance from the central star measured 
on the disc midplane. The initial temperature profile 
corresponds to a flat disc that reprocesses stellar radia¬ 
tion a nd it is steeper when compared with observations 
(e.g. lAndrews et al.l (l2009f ) find a temperature profile 
R~ q ; q ~ 0.4 — 0.74, for discs in Ophiuchus). However 
the actual temperature profile attained once the disc is 
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allowed to evolve is less steep due to additional heating 
sources (viscous and accretion heating). 

We let the disc to relax for 3 kyr 3 outer orbital 
periods) and we then embed in it a protoplanet of ini¬ 
tial mass M p ^ = 1 Mj at distance R p y = 50 AU from 
the central star and let the disc-planet system to evolve. 
The protoplanet’s initial velocity is set the same as the 
velocity of the local gas, i.e. Keplerian (including the 
contribution from the disc mass within the protoplanet’s 
orbit). We assume an initial circular orbit (eccentricity 
ei = 0). The protoplanet is allowed to accrete gas from 
the disc, if this gas approaches within 0.1 AU from the 
protoplanet and is bound to it. This distance is always 
much smaller than the Hill radius of the protoplanet, de¬ 
fined as the region where the protoplanet’s gravity dom¬ 
inates over the gravity of the central star. Therefore, the 
region around the protoplanet is resolved appropriately. 


2.2. Hydrodynamics & radiative transfer 

We use the SPH code seren ( Hubber et al. 2011 ) to 
treat the gas thermodynamics. The code invokes an octal 
tree (to compute gravity and find neighbours), multiple 
particle timesteps, and a 2 nd -order Runge-Kutta integra¬ 
tion scheme. The disc self-gravity is therefore included 
in the simulations. The code uses time-dependent arti¬ 
ficial viscosity with parameters <a m i n = 0.1, <a max = 1 
and ft = 2<a, so as to reduce artificial shear viscosity. 
The chemical and radiative processes that regulate the 
gas temperature are trea ted with t he approximation o f 
iStamatellos et al.l (j2007l ) (see also iForgan et al.l 120091 ). 
We adopt opacity tables appropriate for protostellar discs 
et al. 2M3). 


2.3. Star and protoplanet representation. 

The central star and the protoplanet are represented 
by sink particles that interact with the rest of the com¬ 
putational domain only through their gravity (and lu¬ 
minosity, when irradiation from the protoplanet is taken 
into account). The sink radius of the central star is set to 
A^sink,* = 0.2 AU, and the sink radius of the protoplanet 
is set to R S ink,p = 0.1 AU. This value is always much 
smaller (by a at least a factor of ~ 25) than the Hill 
radius of the protoplanet, defined as the region where 
the gravity of the protoplanet dominates over the grav¬ 
ity of the star, i.e. R S ink, P < Ru = R(M p /3M i ,) 1 ^ 3 . The 
Hill radius increases as the protoplanet accretes material 
from the disc or it decreases as the protoplanet moves 
closer to the central star. Gas particles accrete onto a 
star or protoplanet sink when they are within the sink 
radius and bound to the sink. 


2.4. Radiative feedback from the star and the 
protoplanet. 

The radiation feedback from the star and the pro¬ 
toplanet is taken into account by invoking a pseudo¬ 
ambient radiation field with temperature T A (r) that is a 
functio n of the position relative t o the star and the proto- 
plane t (jStamatellos et a hi 2007, 2011; Sta matellos et all 
l2012f ). This temperature effectively sets the minimum 
temperature that the gas can attain when cooled radia- 
tively. The contribution to T A (r) from the central star is 


set to 

/ r \ —3/4 

T A *(r) = 250K J +10K, (3) 


where R is the distance from the star measured on the 
disc midplane. The contribution to T A (r) from the pro¬ 
toplanet is 


2^planet^ r ) _ f 


L r , 


!6 7T<t sb |r-r p 



( 4 ) 


where L p and r p , are the luminosity and position of the 
protoplanet, respectively. This luminosity is given by 


r ,GM v Mp 
l p — I o 

JVacc 


( 5 ) 


where M v is the mass of the protoplanet, M p is the accre¬ 
tion rate on to it, and R acc the accretion radius. / = 0.75 
is the fraction of the accretion energy that is radiated 
away at the surface of the protoplanet, r ather than be¬ 
i ng e xpended driving jets and/or winds (jMachida et al.l 
2006). As we assume that this planet has formed by 
gravitational instabilities in the disc, the accretion hap¬ 
pens onto the second hydrostatic core. The radius of 
the sec ond core is uncertain; it is estimated to be y 1 — 
20 ([M asunag a fe Inutsukall2000l : iTomida et al.ll20l3l: 
IVavtet et al.ll2013[ ). Here we shall assume i7 acc = 1R 0 . 
The total pseudo-ambient temperature is 

r A 4 (r) = [2»] 4 + [T? lanet (r)] 4 . (6) 

We note that the radiative feedback from the star is fixed 
whereas the radiative feedback from the protoplanet is 
variable with time and depends on the accretion of gas 
onto it. 


3. MIGRATION OF GAS GIANT PLANETS IN 
GRAVITATIONALLY UNSTABLE DISCS 

We investigate the planet-disc interactions and the evo¬ 
lution of the properties of the protoplanet in two cases, 
with and without radiative feedback from the proto¬ 
planet (Figure [1]). The initial evolution of the proto¬ 
planet orbital parameters and mass shows similar pat¬ 
tern in both cases (Figure [2]). The protoplanet ini¬ 
tially migrates inwards fast, with a migration timescale 
of ~ 10 4 years, in agreement with pre vious studies 
(jBaruteau et al.l l20Tll: iMichael et all 1201 lh . As the pro- 
toplanet migrates inwards, its mass increases signifi¬ 
cantly by accreting gas from the disc. Eventually the 
protoplanet is massive enough to be able to open up a 
gap and the migration either slows down (when the ra¬ 
diative feedback from the protoplanet is taken into ac¬ 
count; migration timescale ^ 10 5 years) or even changes 
to a slight outward migration (for the case without ra¬ 
diative feedback from the protoplanet; outward migra¬ 
tion timescale ~ 10 5 years). After the gap opens up, the 
protoplanet continues to grow in mass but rather slowly, 
mainly by accreting material from its circumplanetary 
disc, a disc formed within the Hill sphere of the proto¬ 
planet. The time required for opening up a gap is dif¬ 
ferent between the two cases that we present here. In 
the radiative feedback case the radiation emitted from 
the protoplanet heats the disc, the opening up of the 
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gap is more difficult and takes twice as much time to 
happen (4 x 10 4 yr in comparison with 2 x 10 4 yr for the 
non-radiative protoplanet). The gap not only takes more 
time to open but it is also shallower and less wide (Fig¬ 
ure GBO, as the protoplanet moves closer to the central 
star (where both the disc scale height and Hill radius of 
the protoplanet are smaller). The protoplanet grows in 
mass slower, and the inward migration continues. The 
protoplanet in this case migrates closer to the central 
star (Figure [2^) but its mass growth is considerably sup¬ 
pressed (Figure ( 2 ) 3 ). At the same time the orbit of the 
protoplanet is circularized (Figure^). 

The final outcome of disc-planet interactions is dis¬ 
tinctly different for the two cases presented here (Fig¬ 
ure [T|). In the case without radiative feedback the pro¬ 
toplanet eventually becomes a brown dwarf with a mass 
of 28 Mj, on a eccentric orbit (e ~ 0.17), and a semi¬ 
major axis of 50 AU. When irradiation from the ac¬ 
creting protoplanet is included in the model, the pro¬ 
toplanet reaches a mass of only 14 Mj, i .e. around the 
deuterium burning limit and within the planetary-mass 
regime (< 11 — 16.3 Mi dSoiegel et al.ll2011[) , while it re¬ 
mains on a relatively wide (semi-major axis of 15 AU), 
circular orbit, avoiding excessive inward migration. 

The role of the radiative feedback from the protoplanet 
is critical. Initially, irradiation from the protoplanet 
heats the disc making the opening up of the gap more 
difficult, while thermal pressure delays large gas accre¬ 
tion onto the protoplanet. Therefore, the protoplanet 
continues its inward migration without accreting exces¬ 
sively. Its migration is eventually slowed down once the 
gap is opened up but the protoplanet still continues to 
heat and stabilize the disc. The effective viscosity due to 
the disc self-gravity is small and the inward flow of gas is 
relatively slow. Most of the gas within the protoplanet 
orbital radius is accreted onto the central star creating a 
inner hole in the disc. On the other hand when radiative 
feedback from the protoplanet is not taken into account, 
the disc remains unstable; the Toomre parameter Q at 
the edges of the gap is below 1 (see Figure [3b). There¬ 
fore, the gas flows rather fast towards the outer edge of 
the gap, onto the circumplanetary disc and eventually 
onto the protoplanet, as the effective viscosity due to 
the disc self-gravity is large. The gravitationally unsta¬ 
ble gap edges may also be re sponsible for the outwa rd 
migration of the protoplanet ([Lin fe Papaloizoull2012[ ). 

These results contradict previous findings in which in¬ 
ward migration is fast and continues towards the central 
star as the protoplanet is not able to op en up a gap 
([Barnteau et all l2011t iMichael et al.l [2QTlh . questioning 
whether the survival of giant planets formed early on 
during the lifetime of a disc by fragmentation is possible. 
The critical difference of this work with previous models 
is that the protoplanet is allowed to grow in mass and 
therefore becomes massive enough to be able to open up 
a gap and its migration is slowed down. 

4. NUMERICAL RESOLUTION AND COMPARISON WITH 
PREVIOUS STUDIES. 

To test our numerical code we performed simulations 
of disc-planet interactions in a system where a Jupiter- 
mass planet is embedded in a low-mass (0.005 M 0 ) disc 
at a circular orbit at a distance of 5.2 AU from an 1-M 0 
star. This problem has been studied both numerically, 



0.5 1 1.5 2 

log [column density (g cm' 3 )] 

Fig. 1.— Column density plots showing gap opening induced by 
a Jupiter-mass protoplanet in an 0.1 -Mq disc, in two simulations: 
without (top) and with (bottom) radiative feedback from the pro¬ 
toplanet. The star (in the centre) and the protoplanet are depicted 
by thick white dots. The protoplanet in the first simulation opens 
up a deep, wide gap and grows to become a brown dwarf, migrat¬ 
ing initially inwards and subsequently outwards. In the simulation 
with radiative feedback the gap is shallow and narrow. The proto¬ 
planet migrates inwards but its mass growth is suppressed, so that 
it becomes a wide-orbit planet. 

with grid-based flBate et al.l 120031) and particle -based 
(lAvhffe fc Bat el I2009D codes, and analytically ([Wardl 
1997). We performed simulati ons matching the i nitial 
conditions of previous authors (jAvliffe fe B ate 2009]) us¬ 
ing 10 6 and 2 x 10 6 SPH particles. The Hill radius of the 
planet is resolved adequately: the radius of the planet 
sink is set to 0.1 Ru (where Rn is the Hill radius of the 
planet), while the number of SPH particles that is used 
ensures that the smoothing length, which defines the spa¬ 
tial resolution in SPH simulations, at the planet’s initial 
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Fig. 2. — The evolution of (a) the semi-major axis a, (b) the 
mass M p , and (c) the eccentricity e of the protoplanet. The pro- 
toplanet initially migrates fast inwards, but it grows in mass and it 
is able to open up a gap in the disc. The migration then is essen¬ 
tially halted. In the non-radiative case the protoplanet becomes 
a brown dwarf, whereas in the case with radiative feedback mass 
growth is suppressed and the protoplanet’s mass remains within 
the planetary-mass regime (marked by the horizontal dashed lines 
in (b)). Radiative feedback keeps the protoplanet on an circular 
orbit (e « 0). 


position is 0.17i?H and 0.13i?H for the simulations with 
10 6 and 2 x 10 6 particles, respectively. We find that the 
planet migrates inwards as expected with a migration 
timescale of (1.7 — 2.6) x 10 4 yr when radiative transfer 
is inc lude d (co mpared to 1.5 x 10 4 yr in the literature; 
lAvliffe fe Bate! (j2009D ). and (4 — 7) x 10 4 yr for the local 
isothermal case, i.e. where the disc temperature depends 
only on the distance from the central star (compared to 
(9 — 11) x 10 4 yr and (6 — 12) x IQ 4 yr in the literature; 
lAvliffe Sz Batel ([2009[ ): lBate et al.l (j2003D ). Therefore our 
calculations are in very good agreement with previous es¬ 
timates of migration timescales. 

In the case of the massive, wide-orbit planets that we 
examine here, the planet’s Hill radius is much larger and 
therefore easier to resolve. For example, when the Jo¬ 
vian protoplanet is at 50 AU away from the central star 
its Hill radius is Rn ~ 3.5 AU. We use 10 6 SPH par¬ 
ticles so that the smoothing length at the initial po¬ 
sition of the protoplanet is 0.1 AU, i.e. just 0.03/^h* 
This value remains much smaller than the Hill radius 
as this changes while the planet accretes mass and mi¬ 
grates within the disc (Figure [4]). Another validation 
for our computational method comes from the fact that 
the migration timescale that we obtain for the initial 
stage of planet migration (before a gap is opened up) is 
10 4 yr, which in excellent agreement with pre vious stud¬ 
ies (jBaruteau et alJl201 It IMichael et alJl20lTh using com¬ 


without planet radiative feedback _ (a) 

with planet radiative feedback _ 
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Fig. 3. — (a) Azimuthally averaged surface density of the disc in 
the protoplanet neighbourhood (at 5 kyr ), for the two cases ex¬ 
amined here (without/with irradiation from the protoplanet). The 
radial cells where the surface density is averaged are centred around 
the star, but the distance in the graph is given with respect to the 
cell of the protoplanet (to make comparison easier). The gap is 
shallower and narrower when radiative feedback from the proto¬ 
planet is taken into account, (b) As above, but for the azimuthally 
averaged Toomre parameter Q(R) = c(R) Q(R)/ir G E(R), where 
c is the midplane isothermal sound speed, Q is the the angular 
velocity, E is the surface density, and R is the distance from the 
star. Radiative feedback from the protoplanet stabilises the disc 

(Q > !)• 

plete ly different types of codes: FARGO ([Baruteau et all 
120111 ) and chymera (iMichael et al.ll201lD . 

5. DISCUSSION 

This work demonstrates the combined effect of gas ac¬ 
cretion onto a protoplanet and the associated accretion- 
powered radiative feedback from the protoplanet, on 
regulating its migration and mass growth. Previous 
studies have shown that radiative heating from an 
Earth-like protoplanet forming in a low-mass disc may 
i n fact stop its inward mig ration or even reverse it 
(jBenftez-Llambav et al.l l2015l ). It has also been sug¬ 
gested that radiative feedback may delay the contraction 
of the protoplanetary core so that it i s tidally destroyed 
as it moves closer to the central star (Navaksh in fe Chal 
l2013[ ). The amount of energy radiated from the proto¬ 
planet depends on where the accretion happens, i.e. the 
radius of the protoplanet, which it is uncertain. It has 
been suggested that this energy that is released heats the 
region around the protoplanet, providin g a way to i nfer 
the p resence of a protoplanet in a disc (Montes inos et al.l 
l20l5l) . 

Here we show that accretion and radiative feedback 
work in opposing ways: gas accretion increases the mass 
of the protoplanet and supports the creation of a gap in 
the disc, whereas radiative feedback, which is the result 
of gas accretion, works to suppress accretion onto the 
protoplanet and inhibit/delay gap opening by increasing 
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Fig. 4. — The evolution of the Hill radius, Ru, of the protoplanet 
with time, compared with the smoothing length (0.52 AU), which 
defines the spatial resolution in SPH simulations, at the planet’s 
initial position, and the protoplanet sink radius (0.1 AU). The Hill 
sphere of the protoplanet is resolved adequately. 

the disc temperature. It heats and stabilises the disc, 
suppressing fast gas accretion. Higher gas accretion in¬ 
duces higher radiative feedback that in turn opposes fur¬ 
ther gas accretion. Therefore, accretion and feedback 
are self-regulated and the final outcome of their com¬ 
bined effect is an inward migration of the protoplanet 
with a more moderate increase of its mass than in the 
case without radiative feedback. Eventually the proto¬ 
planet is able to open up a gap, the migration is effec¬ 
tively halted, and the protoplanet survives on a wide 
orbit. 

The final mass of the protoplanet in the simulation 
with radiative feedback is around the deuterium burning 
limit. This simulation covers only a period of 20 kyr, thus 
if gas accretion continues, the mass of the protoplanet 
could further increase. However, radiative feeback may 
delay gas accretion until other disc dispersal mechanisms 
(e.g. photoevaporation) come into play. 


The ramifications of these processes are important for 
the formation of giant planets on wide orbits. Only 
a few planets of this type have been observed so far, 
but this is likely due to observational biases as these 
are faint objects next to bright stars and both high 
sensitivity and high angular resolution observations are 
essential. However, more wide-orbit giant planets are 
bound to be discovered with new surveys using spe¬ 
cialized observing tec hniques (Gemini Pl anet Imager 
([Macintosh et al.l 120141 ). S PHERE/VLT (iBeuzit et all 
l2008f h HiCIAO/SUBARU (jSuzuki et al.l l2009h ) . There- 
fore, it is critical to explain their formation mechanism 
and determine their connection with the formation of 
planetary systems and the formation of our own solar 
system. 

Here we have showed that if massive protoplanets form 
by gravitational fragmentation of discs, gas accretion and 
radiative feedback regulate their mass growth and inward 
migration, so that they are able to avoid excessive mass 
growth and rapid inward migration; their mass remains 
around the deuterium burning limit for at least the ini¬ 
tial stages of their evolution, they avoid disruption and 
they survive on almost circular wide orbits around Sun¬ 
like stars. Therefore, it giant planets like the ones in 
the HR8799 system may have formed by disc fragmen¬ 
tation. Planet formation by gravitational fragmentation 
of protostellar discs may only happen during the early 
stages in a disc’s lifetime, on a timescale of a few thou¬ 
sand years, while discs are still massive enough to be 
prone to gravitational instabilities. The presence of a 
massive giant planet within the disc plays an important 
role in the disc’s evolution and affects the dynamics of 
dust particles as they coagulate to form larger bodies and 
eventually planets within a few million years. Therefore, 
the formation of massive gas giant planets at an early 
stage and their survival on wide orbits sets the disc en¬ 
vironment in which subsequent formation of terrestrial 
planets, like Earth, happens. 
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and from a Royal Society-Daiwa Foundation Interna¬ 
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